Detectability of Gravitational Waves from Phase Transitions 
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Gravitational waves potentially represent our only direct probe of the universe when it was less 
than one second old. In particular, first-order phase transitions in the early universe can generate 
a stochastic background of gravitational waves which may be detectable today. We briefly summa- 
rize the physical sources of gravitational radiation from phase transitions and present semi-analytic 
expressions for the resulting gravitational wave spectra from three distinct realistic sources: bubble 
collisions, turbulent plasma motions, and inverse-cascade helical magnetohydrodynamic turbulence. 
Using phenomenological parameters to describe phase transition properties, we determine the region 
of parameter space for which gravitational waves can be detected by the proposed Laser Interferom- 
eter Space Antenna. The electroweak phase transition is detectable for a wide range of parameters. 

PACS numbers: 98.70.Vc, 98.80.-k 
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I. INTRODUCTION 



Since gravitational waves propagate freely through the 
universe after being generated, their detection provides 
a powerful test of the very early universe. Various mech- 
anisms that generate such gravitational waves have been 
discussed: quantum fluctuations during [l| and shortly 
after inflation Q ; bubble wall motion and collisions dur- 
ing phase transitions @, El IS IS 0i cosmic strings and 
defects, including primordial black holes H; cosmolog- 
ical magnetic fields [S US El, and plasma turbulence 
[H, UH [l5j] . This paper focuses on gravitational waves 
generated during phase transitions and study the possi- 
bility of direct detection by the planned Laser Interfer- 
ometer Space Antenna (LISA), 

Near-future data from the Large Hadron Collider will 
for the first time have the ability to probe in detail 
physics at the electroweak energy scale. This physics de- 
termines the nature of the electroweak phase transition 
in the early universe, when the primordial plasma went 
from an electroweak-symmetric state to a broken state 
with distinct electromagnetic and weak interactions. The 
phase transition took place when the primordial plasma 
had a temperature on the order of 1 TeV. Intriguingly, 
the Hubble frequency f/» at this epoch, redshifted to to- 
day, falls in the lower end of the detection range for LISA, 
which ranges from 1CP 4 to 1CP 1 Hz [la ]. If gravitational 
waves are generated during the electroweak phase transi- 
tion, their characteristic frequency will be related to this 
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Hubble frequency: remarkably, gravitational waves may 
provide us with an alternate route to probing electroweak 
physics. 

Having gravitational waves at detectable frequencies 
is not sufficient for detecting the electroweak phase tran- 
sitions: the source of gravitational waves must be suf- 
ficiently strong to produce radiation with a detectably 
large amplitude today [13, ES ES H3] • If the phase tran- 
sition is first order, the latent heat of the phase transition 
partly is transferred into kinetic energy of the walls of ex- 
panding bubbles of the broken phase. If these expanding 
bubbles contain a large enough amount of energy, they 
will produce a gravitational wave background with a de- 
tectably large amplitude today [3]. Models of the elec- 
troweak phase transition based on standard-model parti- 
cle physics do not produce observable gravitational wave 
signals since they are not first order for allowed values 
of the Higgs mass [5l[, but common extensions of the 
standard model, including supersymmetry and extra di- 
mensions, can produce much stronger phase transitions 

Bubbles of broken phase in a first-order phase transi- 
tion expand and percolate to convert the entire universe 
to the low-energy phase. The kinetic energy in the bubble 
walls eventually thermalizes, but prior to that the bub- 
bles act to stir the primordial plasma, plausibly gener- 
ating Kolmogoroff turbulence cascading from the turbu- 
lence scale to a much smaller scale where the kinetic en- 
ergy turns into thermal energy via viscous heating. The 
turbulent eddy motions can also be a potent source of 
gra vitational waves if the energy input is large enough 
[l2j . If a small magnetic seed field with nonzero helic- 
ity is generated at the phase transition [23|, [24j] , magne- 
tohydrodynamic effects can generate an inverse cascade, 
transferring energy to scales substantially larger than the 
stirring scale [HI and resulting in a detectably large sig- 
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nal [25j], and additionally can give the radiation a non- 
zero circular polarization [l4j] which is also directly de- 
tectable [26[. So gravitational wave generation from an 
early universe phase transition can be decomposed into 
three distinct sources: expanding bubbles of the broken 
phase, hydrodynamic turbulence stirred by the colliding 
bubbles, and an inverse cascade due to the amplification 
of seed magnetic fields. 

Previous papers have computed the gravitational wave 
spectra for each of these sources individually. Here we 
consider their combined spectrum. The detailed shape 
and amplitude of the spectrum depends on fundamental 
properties of the phase transition: its energy scale, bub- 
ble nucleation rate, latent heat, efficiency of converting 
latent heat into plasma motions, and mean helicity of 
seed magnetic fields. We summarize how the total power 
spectrum depends on these parameters, then evaluate the 
regions of parameter space for which the relic gravita- 
tional radiation from the phase transition is detectable 
with LISA, displaying the results as contour plots in the 
parameter space. While the bubble collisions produce de- 
tectable gravitational radiation for a very limited range 
of parameters, adding the radiation from turbulence and 
an MHD inverse cascade greatly widens the range of de- 
tectable phase transitions. 

In the following Section, we give a brief overview of the 
physics of phase transitions and the corresponding phe- 
nomenological parameters. In Sec. Ill, we express the 
gravitational wave spectra in terms of these phenomeno- 
logical parameters, and show how the total gravity wave 
spectrum varies with the parameters. Sec. IV compares 
these spectra with the projected LISA sensitivity curve 
for stochastic backgrounds, displaying detectability re- 
gions in the space of phase transition parameters. 



II. MODELING GRAVITATIONAL RADIATION 



phase transition lasts for a finite duration, we can con- 
sider the total radiation field to be an incoherent sum 
of radiation from many sources, each with a size corre- 
sponding to the light travel distance during the phase 
transition. With this approximation, we can apply the 
far-field approximation for each individual source region, 
x 3> d where d is the size of the source region; in this 
region the gravitational waves are the only metric per- 
turbations [281] . allowing replacement of |x — x'| by |x| in 
Eq. CQ. 

Generated gravitational waves propagate freely 
through the expanding universe until today. Their 
wavelengths simply scale with the scale factor a of the 
universe, while their total energy density evolves like a -4 
and their amplitude decays like a -1 . Following Ref. [27], 
we use the gravitational wave spectral energy density 
Pgw{u) per logarithmic angular frequency w, defined 
as the root-mean-square average of pcwfejUj) over 
spatial positions, and form the spectral energy density 
parameter Q,qw{u) — Pgw(u)/ p c with the critical 
energy density p c = 3Hq/8ttG. Changing to linear 
frequency / = lu/2tt, a characteristic strain amplitude is 
conventionally defined as 



h c (f) = 1.3 x 10 



-18 



1Hz 
~T 



[^fi GW (/)] 1/2 



(2) 



where h$ is the current Hubble parameter Hq in units 
of 100 kmsec _1 Mpc _1 . We will present our results in 
terms of the gravitational wave amplitude h c vs. linear 
frequency / measured today. 

The temperature at which the phase transition occurs 
determines the corresponding Hubble frequency, which is 
given today by [IH 



/ H = 1.6 x 10 _5 Hz 



T* 



lOOGeV/ V 100/ 
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(3) 



A. General Remarks on Gravitational Waves 
Generation 

Gravitational waves generated during phases transi- 
tions by stochastic sources have an energy density just 
after the phase transition of [27] ] 



PGw(x) 

G_ 

~ 2tt 



1 
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x — x x 



(1) 



where time variables with primes represent the light- 
cone combination t' = t — |x — x'|, i and j are spa- 
tial indices (repeated indices are summed), the source 
Sij(x,t) = Tij(x,t) — 5ijTj:(x,t)/3 is the traceless part 
of the stress-energy tensor Ty, G is the gravitational 
constant, and we use natural units with h = c = 1. 
This expression assumes that the duration of the phase 
transition is short enough to neglect the cosmological ex- 
pansion during gravitational wave production. Since the 



where is the number of relativistic degrees of freedom 
at the temperature T*; for standard model degrees of 
freedom, ~ 106.75 as T — > oo. To obtain this expres- 
sion, begin with the Hubble parameter at temperature 
T*, given by 



m = 



8tt 6 T? 
90M|, 



(4) 



then rescale H* to its present value with the ratio be- 
tween the scale factors a* corresponding to the Universe 
temperature T* and ao today: 



(5) 



Note that an extra factor of approximately 3.6 occurs 
for /h in Ref. [1§] (Eq. 10) due to approximating H% ~ 
Gptot- 

It can be shown from general physical principles that 
since the gravitational wave generation process is causal, 



3 



the peak frequency must be greater than or equal to the 
Hubble frequency, / pea k > /h |. Independent of the 
nature of the source, the gravitational wave spectrum is 
the same at low frequencies, flaw if) cx (/// pea k) 3 (for 
/ « /peak) @, HO, M, M, M HI The gravitational 
spectrum behavior at higher frequencies is determined 
by the specific source features. 



B. Phase Transition Model Parameters 



A first-order phase transition is generically described 
by several parameters: (i) a = /0 V ac/pthcrmai, the ratio of 
the vacuum energy associated with the phase transition 
to the thermal density of the Universe at the time (which 
characterizes the strength of the phase transition); (ii) k, 
an efficiency factor which gives the fraction of the avail- 
able vacuum energy which goes into the kinetic energy of 
the expanding bubble walls, as opposed to thermal en- 
ergy; (in) which sets the characteristic time scale for 
the phase transition; (iv) Vb, the velocity of the expand- 
ing bubble walls, which set the characteristic length scale 
of the phase transition; (v) T*, the temperature at which 
the phase transition occurs. In any first-order phase tran- 
sition, the characteristic bubble nucleation rate per unit 
volume is generically [3l| 



r = r oe 



lit 



(6) 



which uses the constant and linear terms of a Taylor ex- 
pansion of the bubble nucleation action. General consid- 
erations confirmed by numerical calculations show that 
the largest bubbles reach a size of order j3~ 1 Vb by the end 
of the phase transition [6( , where Vb is the bubble expan- 
sion velocity, assuming the bubbles remain spherical as 
they expand. In general, j3 is expected to be of the order 
41n(rn P i/T)i? ~ 100H for a Hubble rate H 

The fundamental symmetry breaking mechanism 
which drives the phase transition determines some effec- 
tive potential for bubble nucleation [32| . The difference 
in energy density between the two phases and the bubble 
nucleation rate are both determined by this mechanism. 
Thus the parameters T* , /?, and a are all determined di- 
rectly by the underlying physics. On the other hand, the 
bubble velocity i>b and the fraction of kinetic energy into 
the bubbles k depend on the detailed microphysics in- 
volved in the bubble propagation through the relativistic 
plasma and are not determined from general properties 
of the effective potential. In general, the larger the vac- 
uum energy density driving phase transition, the higher 
the bubble wall velocity Vb 0, Hl| . 



III. SOURCES OF GRAVITATIONAL 
RADIATION 

A. Bubble collisions 

Bubbles walls can propagate outwards through a rela- 
tivis tic p lasma via two modes, detonation and deflagra- 
tion [33|. In the case of detonation, the bubble walls are 
thin compared to the bubble radius with velocity is [H| 



v b {a) = 



l/v / 3 + (a 2 + 2a/3) 1 / 2 
1 + a 



and an approximate form of n is given 



k(o) 



1 



1 + Aa 



Aa 
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(7) 



(8) 



with A — 0.715. If the bubbles propagate as a deflagra- 
tion front, the walls are thick and have lower energy den- 
sity. However, it has been argued [34| that for relativistic 
plasma, instabilities in the bubble shape will accelerate 
the bubble walls and the hydrodynamic deflagration ex- 
pansion mode is unstable to becoming a detonation. For 
this reason, Ref. [1] assumed that Eqs. (0) and |J5]) hold, 
and we will also make that assumption here. (For unusual 
cases with very strong detonation and a > 1, the follow- 
ing formulas must be corrected.) Both modes of bubble 
propagation recently have been re-addressed in Ref. [3(| , 
with results slightly different from those of Ref. [H . 

The peak frequency of the gravitational wave spectrum 
is of fundamental importance for issues of detectability. 
The bubble nucleation time scale is always determined 
by /3 _1 , while the largest bubble size depends also on the 
bubble wall velocity, Iq = Vbfl^ 1 - What scale determines 
the peak frequency of the gravitational waves? When 
the bubble wall velocity Vb approaches the speed of light, 
Iq 1 ~ /?, and both characteristic frequency scales are the 
same. In Refs. [1, IH @, where the bubble wall velocity is 
assumed to be the speed of light, the peak frequency is 
given by 0, MM 



/peak = 0-3 I ] / H 



/ peak 



(9) 



Note the angular peak frequency at the moment of emis- 
sion Upc &kir is set by 0.3k Vb ~ 2(3 (where fc = 2-k/Iq). 
To understand the peak frequency of the radiation from 
bubble walls with lower velocity, it is useful to exploit the 
direct analogy with sound wave generation. The typical 
source velocity compared to the velocity of the induced 
waves determines the characteristic peak frequency [35j. 
If the characteristic source velocity, given by the ratio be- 
tween the source's characteristic length and time scales, 
is less than the group velocity of the induced waves, then 
the frequency is determined by the source time scale; 
in the opposite case, the characteristic length scale de- 
termines the peak frequency. Applied to gravitational 
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waves with group velocity v = 1, the source time scale 
always determines the peak radiation frequency, since by 
causality its characteristic velocity can never exceed the 
speed of light. 

For the realistic values ft/H* — 100 and = 100, the 
peak frequency from bubble collisions /p Cak — 5 x 10~ 4 Hz 
[HI is close to the peak of LISA sensitivity at 2mHz if 
T* ~ lOOGeV. The corresponding peak amplitude of the 
gravitational radiation today is [4| 



M/pcak) = 1-8 X 1(T 16 A 



n 



/100Gev\ 



a + 1 

5 







100 

9* 



0.24 



(10) 



The spectrum shape at high frequencies (/ ^> /h) is 
h c (f) « f~ 2 ver y cl° se to the scaling h c (f) oc / _19 
found numerically in Ref. [ij. 



B. Turbulence 

Once the bubbles expand and percolate, a significant 
amount of their kinetic energy will be converted to tur- 
bulent bulk motions of the primordial plasma Q . As the 
phase transition ends, far more small bubbles have been 
nucleated than large ones, but the energy density in the 
large ones dominates the total energy density [3l[ . There- 
fore, a reasonable approximation is that the turbulence 
energy is injected on a stirring scale lo corresponding to 
the size of the largest bubbles, and that that the stir- 
ring will last for roughly T st i r = comparable to the 
duration of the phase transition. We make the conser- 
vative assumption that the duration of the turbulence is 
also r s ti r . The turbulence can actually last substantially 
longer than this; see Ref. [12] for a detailed discussion. 
After generation, the turbulence kinetic energy cascades 
from larger to smaller scales, stopping at a viscous damp- 
ing scale Id where the turbulence energy is dissipated into 
heat. We define the parameter 7 = IqH^ ~ Vb{H*/j3) 
which determines the number (iVeddy — 7 -3 ) of turbu- 
lent eddies within the Hubble radius. In terms of the 
parameter 7, Eq. © reads = 0.34/ H V7- 

We associate a characteristic velocity perturbation vq 
with the fluid velocity of the largest eddies at the stirring 
scale lo. This velocity determines the turbulent motions' 
Mach number M ~ vq- It is easy to estimate vo through 
relating the kinetic energy density of plasma w«g/(l — Vq) 
(where w = p + p is enthalpy of plasma) and the vacuum 
energy p vac using w = iap vac /3. This leads to [l8T ] 



4 + 3Ka' 



(11) 



Accounting for Eq. |(5J), vq can be expressed uniquely in 
terms of a and in the case of realistic phase transitions, 
when a < 1, Vq ~ a 3 ^ 4 (3/2) 1 / 4 /3, and of course is less 
than the bubble wall expansion speed Vb- For kq ~ 1 



corresponding to a strongly first order phase transition, 
v ^ 0.65. 

Figure 1 of Ref. [l5[ shows that it is safe to adopt the 
aero-acoustic limit k — > in the gravitational wave gen- 
eration formalism for turbulent sources, provided only 
that the Mach number M < 1. Employing this simplify- 
ing assumption and adopting the results of Ref. [15j , the 
peak frequency of gravitational waves is given by 



J peak 



I /h 1 T T 



3 I -V 



Vb J 



peak - 



(12) 



For typical parameters (3 = lOOii*, T* = 100 GeV, g* = 
100, and a = 0.5 (which correspond to ft = 0.36, Vb = 
0.89 and vo = 0.34), the peak frequency is 0.4 mHz, 
within the LISA range. 

At frequencies high compared to the characteristic 
damping frequency, / > f D = /J cak i? 1/2 with R is 
the turbulence Reynolds number, the gravitational wave 
strain spectrum h c {f) ~ (///d) 1/2 exp(-(f / f D ) 2 ) pos- 
sesses a sharp exponential cutoff. Near the peak fre- 
quency, for fp Cak < f < Id, the spectrum scales like 

hc(f) pi (J7JjL k r 13/4 El- The peak frequency of the 
radiation spectrum from bubble collisions for the param- 
eters considered here is lower than the peak frequency 
due to turbulence, fp^ k < /J oa k, and the bubble col- 
lision spectrum falls only as {f / f^cak) 2 at frequencies 
above the peak, but even at frequencies substantially 
greater than /^ ak , the total gravitational wave spectrum 
is not dominated by the bubble collision source, due to 
the lower amplitude of gravitational waves induced by 
collisions. 

The gravitational wave amplitude from turbulence can 
be conveniently estimated analyti cally using the aeroa- 
coustic approximation (see Refs. jl5l . |25| for details of 
numerical computations), giving 



a + V 



'100GeV\ 

[ — ~ 1.9M2 



v 3 h +0.24 J kCfp.Sk) 



(13) 



and for a ~ 0.5 is approximately 3 times larger that the 
peak amplitude coming from bubble collisions. 



C. MHD Inverse Cascade Turbulence 

Several models lead to generation of a cosmolo gica l he- 
lical magnetic field during phase transitions [23|, l24( . In 
such a case, turbulence develops in a magnetized plasma, 
which is qualitatively different than an unmagnetized 
plamsa. In particular, an inverse cascade which gener- 
ates magnetized perturbations on larger scales than the 
stirring scale can occur. To model such turbulence, we 



assume that the hydrodynamic and magnetic Reynolds 
numbers are much greater than unity on scales ~ Iq 
(which is simply the condition for turbulence to develop) ; 
we also assmue that the dynamics of magnetohydrody- 
namic turbulence is dominated by Alfven waves for which 
the magnetic and kinetic energy densities are in approx- 
imate equipartition [36]. We also assume a small initial 
magnetic helicity, C* = Hm(^/[2£m(**)£m~(**)] < 1, 
of a magnetic field with the magnetic-eddy correlation 
length f M (t) = J E M (k,t)k~ 1 dk/£ M (t) and with to- 
tal magnetic energy £m(0 = / E M (k,t)dk and helic- 
ity Ti.M(t) = J H M (k,t)dk. From the stirring scale, 
the direct cascade proceeds via a Kolmogoroff-like model 
[37L HH and the following inverse cascade stage adopts 
the MHD turbulence model of Refs. jjH, Hj. Then it 
can be shown that the induced gravitational wave spec- 
trum always peaks at the Hubble frequency = fa 
[25l | . and compared to the turbulence peak frequency it 
is shifted to lower frequencies by a factor (vb/M){H*/ (}). 

The amplitude is larger compared to hydrodynamic 
turbulence case as [H} 

= (^) 3/4 C, 9/8 / ic (/ p T cak ). (14) 

This increase in amplitude is because the duration time 
of the inverse cascade is longer than the duration of the 
direct cascade. No efficient dissipation mechanisms ex- 
ist at the largest scales, so the cascade will stop at a 
scale £m = 1/ks either when the cascade time scale r cas 
reaches the expansion time scale H~ l or when the char- 
acteristic length scale £m(£) reaches the Hubble length 

H^ 1 . These conditions are 1 ^1%/vqIq < H^ 1 or 
lg = 27r/fcs < H^ 1 (the cascade time is scale depen- 
dent and maximal at k = ks)- Defining 7 < 1, it is 
easy to see that the first condition is fulfilled first, and 

consequently k Q /ks < (vq/j) 1 ^ 2 C*^ • To have an inverse 
cascade requires k^/kg > 1, leading to a constraint on 
1/2 

initial helicity 7 < MC,J . In addition, these gravita- 
tional waves will be circularly polarized since they are 
induced by parity- violating stochastic sources [T3| • If al- 
ternately we use the inverse cascade model of Ref. [4lj |. 
the peak frequency of the gravitational wave spectrum 
remains unchanged while the peak amplitude is doubled 

USUI. 



D. Total Gravitational Wave Spectrum 

The total gravitational wave background from a first- 
order phase transition is a sum of three terms correspond- 
ing to the three distinct sources discussed here. Our re- 
sults are expressed in terms of the effective strain spec- 
trum h c (f) in terms of the phase transition temperature 
T* , the number of effective relativistic degrees of freedom 




1Cf 6 10" 5 1CT 4 10" 3 1Cf 2 

f/Hz 



FIG. 1: The spectrum of gravitational radiation for a first- 
order phase transition with g, — 100, T„ = 100 GeV, a = 0.5 
and j3 = 100-ff*, from bubble collisions (dash line), hydro- 
dynamic turbulence with zero helicity (solid line) and MHD 
turbulence with f» = 0.15 (dash-dotted line). The bold solid 
line corresponds to the 1-year, 5a LISA sensitivity curve 14311 . 
including the confusion noise from white dwarf binaries [441 ] . 
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FIG. 2: The total spectrum of gravitational radiation (in- 
cluding there sources of gravitational radiation bubble colli- 
sions, hydro-turbulence, and MHD turbulence) for = 100, 
T, = 100 GeV, C* = 0.1, (3 = 100#*, and three different val- 
ues of a: a = 1 (solid line), a = 0.5 (dash line) and a = 0.2 
(dotted line), with the LISA sensitivity curve.. 
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FIG. 3: The total spectrum of gravitational radiation for 
g* = 100, a = 0.5, = 100if„, C* = 0.1, and three different 
temperature values, T» = 100 GeV (solid line), T, = ITeV 
(dash line), and T» = lOTeV (dotted line), with the LISA 
sensitivity curve. 
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FIG. 4: The total spectrum of gravitational radiation for 
g, = 100, T» = 100 GeV, a* = 0.5, (3 = 10077*, and four 
different values of f«: = 0, corresponding to hydrody- 
namic turbulence without the inverse cascade effect (solid 
line), C* = 0-02 (dash line), = 0.5 (dotted line), and 
f* = 0.15 (dot-dash line), with the LISA sensitivity curve. 



r?*, the ratio of vacuum energy to thermal energy a, the 
bubble nucleation time scale /?, and the initial magnetic 
helicity parameter £* . All others relevant parameters ap- 
pearing in the analytical approximations, such as k, Vb, 
Iq, and M can be computed in terms of these fundamen- 
tal parameters using simple and well-motivated assump- 
tions. We fix g* = 100, and we do not study the de- 
pendence on this parameter since it is very weak and the 
actual value of <?* will not vary more than a factor of a few 
from this nominal value (g* = 106.75 for standard model 
degrees of freedom at temperatures large compared to the 
W and top quark masses). We choose (* small enough 
to insure that the resulting magnetic turbulence can be 
modeled as an inverse cascade modeling (see Ref. [25j for 
details). 

Figure 1 displays the three components contributing 
to the total gravitational wave spectrum, for a fiducial 
model with g* = 100, a = 0.5, fi/H* = 100, andT* = 100 
GeV. Bubble collisions are shown as the dashed line, hy- 
drodynamic turbulence with no magnetic helicity as the 
solid line, and MHD turbulence with an initial helic- 
ity parameter £* = 0.15. The dark curve is the LISA 
design sensitivity curve for stochastic backgrounds with 
a detection signal-to-noise ratio of 5 over a 1-year in- 
tegration [43)], including a rough estimate of the con- 
fusion limit from white-dwarf binaries [44| . It is clear 
that bubble collisions alone do not produce a significant 
signal (marginally detectable) in LISA for this range of 
phase transitions parameters, agreeing with previous re- 
sults [ril [iH [T9L |30]| . The turbulence peak is observable, 
as well as the separate peak from the MHD inverse cas- 
cade. 

The next set of figures show how the total gravitational 
wave spectrum varies as each of the phase transition pa- 
rameters is changed, holding the others fixed. Figure [2] 
shows the variation with a, Fig. [3] the variation with T*, 
Fig-Sthe variation with £*, and Fig.[5]the variation with 
P/H+. In each case, the spectrum for the fiducial model 
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FIG. 5: The total spectrum of gravitational radiation for g, — 
100, T* = 100 GeV, C» = 0.1, a = 0.5, and three different 
values of /3: (3 = 40//, (solid line), /3 = 10OH, (dash line) and 
(3 — 500//* (dotted line), with the LISA sensitivity curve. 



in Fig. 1 is plotted, along with the spectra for two or 
three other values of the parameter in question. We can 
see qualitative outlines of detectability from these figures. 
For example, Fig.[3]shows that as the temperature of the 
phase transition increases, detection relies increasingly 
on the MHD signal. 

While these plots give a good qualitative sense of 
the spectrum dependence on the parameters, we would 
like to know more precisely which regions of the four- 
dimensional parameter space of phase transitions result 
in a gravitational wave spectrum detectable with LISA. 
The following section gives one useful representation of 
the detectability regions. 



IV. DETECTABLE REGIONS IN PARAMETER 
SPACE 

The remaining figures show various detectability re- 
gions. For the first set of four plots, Figs. [6] to [9j the 
detectable region in the [3/H+-a plane is displayed; each 
figure is for a different value of the phase transition tem- 
perature T*, ranging from 0.1 TeV to 100 TeV. The larger 
displayed region is including the MHD inverse cascade 
from an initial helicity of £ = 0.15, while the smaller re- 
gion shows the case of no inverse cascade, or equivalently 
the limit ( = 0. The other parameter, g„, is held fixed 
and only has a weak effect on the results. We designate 
a power spectrum as detectable if its amplitude comes 
above the LISA 1-year sensitivity curve at signal-to-noise 
ratio of 5 at any frequency; more careful detection con- 
ditions can also be used but will give the same general 
results. 

Large regions of parameter space are detectable. The 
natural expectation for (3 is around lOOiZ* 31], with some 
range around this value for plausible models. For /3/iJ* 
between 10 and 100, the detection threshold in a tends 
to be roughly constant, with lower a being detectable as 
T* increases. If T* > 1 TeV, then phase transitions with 
a > 0.06 are detectable by LISA. 
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FIG. 6: The a vs. /3/H* LISA sensitivity region for a phase 
transition at a temperature T* =0.1 TeV. The regions for 
C» = and = 0.15 coincide at this temperature. A point in 
parameter space is considered detectable if at any frequency 
its value of h c (f) is detectable at a signal-noise ratio of 5 in a 
one- year integration, including the confusion noise from white 
dwarfs from 14311. 




FIG. 7: Same as Fig. El except for T* = 1 TeV. For C» = 0.15, 
a slight additional region along the low-a edge of the = 
region is detectable. 



For temperatures < 1 TeV, the contribution from 
the MHD inverse cascade is not important for LISA de- 
tectability, due to the experiment's low-frequency cut- 
off. As T* moves to higher temperatures, the typical fre- 
quency, set by ii*, becomes larger, and the low- frequency 
peak and tail from the MHD inverse cascade starts to 
come into the detection window. For T* > 10 TeV, the 
detectability range in (3 is substantially expanded due 
to the magnetic source, though mostly at comparatively 
large values of a. The spectrum contribution from hy- 
drodynamic turbulence, due to its somewhat higher fre- 
quency and amplitude than that from bubble collisions 
and much higher frequency than the MHD contribution, 



FIG. 8: Same as Fig.U except for T* = 10 TeV. For C* = 0.15, 
a substantial extra area on the right side of the = region 
is detectable. 




FIG. 9: Same as Fig. except for T, = 100 TeV. 



dominates most of the detectability for the LISA fre- 
quency band. 

Figures flOl and fTTI show the detectability regions in the 
/3-T* plane, for a = 0.1 and a = 0.5, respectively. For 
a = 0.1, a relatively large region is detectable, while for 
a = 0.5, a phase transition at any temperature up to 10 
TeV is detectable for any reasonable value of f3. 

These results will be useful for computing future con- 
straints on various high-energy physics theories which 
give first-order electroweak phase transitions. The funda- 
mental parameters of any theory can be used to construct 
order parameters and vacuum energies for any phase 
transition at relevant energy scales (see, e.g., Refs. [22[ 
for many example). These in turn can be used to com- 
pute the phcnomcnological phase transition parameters 
a, (3, and T* used here, as done in Appendix A of Ref. [|[. 
The parameter likely will require additional astrophys- 
ical modelling in situations where there is no well-defined 
mechanism for generating magnetic fields, but most of 
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FIG. 10: The /3/H* vs. T, LISA sensitivity region for a phase 
transition with vacuum energy a = 0.1. The regions for = 
and C,* = 0.15 coincide at this temperature. A point in 
parameter space is considered detectable if at any frequency 
its value of h c (f) is detectable at a signal-noise ratio of 5 in a 
one- year integration, including the confusion noise from white 
dwarfs from 14311. 




FIG. 11: Same as Fig.[TDl except for a = 0.5. 



the parameter-space detectability is independent of this 
parameter. 

We have made several basic physical assumptions: we 
have assumed that the walls of the bubbles formed in the 
first order phase transition expand as detonations [34| ; 
that the bubbles remain spherical as they expand; that 
the resulting bubble collisions source turbulence on the 
scale of the largest bubbles for a time comparable to the 



duration of the phase transition; that the turbulence can 
be modeled as fully-developed Kolmogoroff turbulence; 
and that initial seed magnetic fields with non-zero helic- 
ity will be amplified in an inverse cascade according to 
a specific MHD turbulence model [39j, |40j • While these 
are all reasonable assumptions and are likely close to the 
realistic physical situation, they are not guaranteed to 
be correct. In particular, the assumption that bubbles 
in a relativistic plasma always expand as detonations de- 
serves further study. The transition from coherent bub- 
ble wall expansion to stochastic turbulent behavior needs 
to be simulated numerically; current computing power is 
sufficient to model this process with enough resolution to 
verify the simple assumptions in this work. The standard 
Kolmogoroff turbulence results apply to non-relativistic 
fluids, and our direct extension to relativistic plasmas, 
while likely qualitatively correct, is not based on either 
direct observation or detailed calculation. Finally, mech- 
anisms for generating seed magnetic fields prior to or dur- 
ing the phase transition are not clearly understood. The 
enticing prospect of a direct probe of electroweak physics 
via the stochastic background of gravitational radiation 
is a strong motivation for further investigation of these 
interesting questions concerning the fundamental physics 
of the early universe. 
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